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Abstract 

We report on the discovery of new Martian Trojans within the Minor Planet 
Center hst of asteroids. Their orbital e volution over lO^yr shows character- 
istic signatures of dynamical longevity (jSchoU et al.l . |2005[ ) while their aver- 
age orbits resemble that of the largest known Martian Trojan, 5261 Eureka. 
The group forms a cluster within the region where the most stable Trojans 
should reside. Based on a combinatorial analysis and a comparison with the 
Jovian Trojan population, we argue that both this feature a nd the appar- 



ent paucity of km-sized Martian Trojans (ITrilling et al. 



to expectations from earlier work ( iTabachnik and Evans 



2006 ) as compared 



19991 ) is not due to 



observational bias but instead a natural end result of the coUisional comminu- 



tion 


Jutzi et al. 


2010 


of a proge 



5r L5 Trojan of Mars. Under the collisional scenario in 
particular, the new Martian Trojans are dynamically young, in agreement 
wi th our age es timate of this "cluster" of < 2 Gyr based on the earlier work 



of IScholl et al.l . This work highlights the Trojan regions of the Terrestrial 
planets as natural laboratories to study processes important for small body 
evolution in the solar system and provides the first direct evidence for an 
orbital cluster of asteroids close to the Earth. 
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1. Introduction 



Trojan asteroids are objects confined by solar and planetary gravity to 
orbit th e Sun 60° ahead (La) or behind (L5) a planet's position along its 



orbit (IMurray and Dermottl . Il999l ) . The y hold a special status among aster- 



oids. D u e to their dynamic al longevity (ILevison et al.l . 119971 : iMarzari et al. 



2003al Jbl: IScholl et al.l . |2005| ). they are thought to be left o yer material from 



the formation and early eyolutiq n of our planetary system ([Shoemaker et al 



Nesvorny and Vokrouhlickyl . l2009l : iLykawka and Hornerl . I2OIOI ) . Mars is the 



19891: iMarzari and Scholill998alJbl:lMarzari et al.l.l2002l:lMorbideUi et al.l . 12005 



only t errestrial planet cur rently known to host stable Trojans (jSchoU et al 



20051 : iDvorak et al.l . l2012l ). As of 21 June, 2012, only three long-lived ob- 



jects appear in the Minor Planet Center (MPC) list of Martian Trojans 



(http:/ /www.minorplanetcenter.net/iau/lists/MarsTrojans.html|: 5261 Eu- 
reka (L5), (101429) 1998 VF31 (L5) and (121514) 1999 UJ7 (L4). Inclu- 
sion in this list requires that they be confirmed as Trojans by numerical 
integrations for at least 10^ yr (Gareth Williar ns, priv. com m.). Several 



other objects are temporary co-orb it als of Mars (IConnors et al.l. l2005l ) , sim- 



ilar to Earth's transient co-orbitals (iNamouni et al. 



Morals and Morbidellil . 120021 : iBrasser et all 120041 : IWaierl . I2OIOI ). The debate 



1999 



Christou. 2000 



of the existenc e and origin of long- lived terrestrial planet Trojans is now com- 
ing to the f ore (Todd et al. . 2012al b) due to expecte d results from the up com- 
ing GAIA (iMignard et al.l . l2007f l and NEOSSSAT (llaurin et all . I2OO8I ) mis- 
sions as well as new deep, wide field surveys of the sky such as the Panoramic 



Surve y Telescope & Rapid Response System (PanSTA RRS: jjedicke et al. 



20071 ) and the Large Synoptic Survey Telescope (LSST: jjones et al.l. 120091) . 



A Trojan of the Earth was discovered rec entlv (IConnors et al.l . 1201 ll ) and 



long-lived Trojans of our planet may exist ( Dvorak et al. . 20121 ). 

The expected n umber of Martian Trojans > 1 km in diameter ha s been es- 

timatedtob e < 50 (iTabachnik and Evansl . ll999l ). Previous studies ( ITabachnik and Evans . 



19991 . l2000al : IScholl et al.l . l2005l ) have shown that the most long-lived Tro- 



an librators of M ars exist in particular regions of orbital element space. 



Scholl et al.l (120051 hereafer referred to as SMT05) mapped out Trojan sta- 



bility as a function of their proper elemen ts, quantities that are constants 
ofmotion under planetary perturbations ( Milanil . 1993 : Beauge and Roig , 
2OOII ). They found that the most stable objects are those with proper ec- 



centricity Bp < 0.15, proper inclination 13° < Ip < 28° and proper libration 
amplitude 0° < D < 150°. Eureka and 1999 UJ7 (e^ ~ 0.05, Ip ~ 22°) 
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occupy the projection of this stable island on the Cp — Ip plane, whereas 1998 
VFsi's higher proper inclination (Jp ~ 30°) places it in a different region, 
separated from the stable island by the z/5 secular resonance. Both 1999 
UJ7 and 1998 VF31 have proper libration amplitudes > 50°, unlike Eureka 
(~ 10°). Further, SMT05 found that long-lived Trojans (ie not the transi- 



tory types investigated by IConnors et al.l ( 120051 )) belong to two groups: the 
"fast diffusers" (FD) that escape from Trojan libration within 10*^ yr and 
the "slow diffusers" (SD) that persist as Trojans of Mars for periods of time 
comparable to the age of the solar system. 

Here we report on a study of new candidate Martian Trojans discovered 
within the Minor Planet Center list of asteroids. The osculating orbital 
elements of these objects are similar to Eureka's, suggesting that they may 
also be long-lived. This work aims to determine, in the first instance, whether 
these are transitory, and if not, whether they belong to the SD or the FD 
group. Additionally, we wish to constrain the nature of their relationship, 
if any, to Eureka and the other known SD Trojans. In the next Section we 
describe the new Trojan candidates. In Sections [3] & IH we establish their 
stability and relation to the other known Trojans. In Section [S] we quantify 
the statistical significance of the observed orbital clustering while in Section [6] 
we argue that it is likely not an artifact of observational bias. In Section [7] 
we compare the Martian and Jovian Trojan populations. In Section 8 we 
place our findings in context and discuss different scenarios for the origin 
and evolution of these objects to the present date. In the last Section we 
summarise our main conclusions and outline the work's implications for small 
body research. 

2. Search for new candidate Martian Trojans 

We searched through the Minor Planet Center databas^ and found 7 
asteroids with similar osculating semimajor axis, eccentricity and inclination 
to Eureka. The specific criteria used were: 1.523 AU < a < 1.5245 AU, 
e < 0.15 and 16° < I < 24°. Table [1] shows these eight objects and, in 
addition, the long-term stable Trojan with significantly higher inclination 
than Eureka, 1998 VF31. We consider the six objects that are not mentioned 
in the MPC list to be new candidate Martian Trojans. 



^Available at http://www.ininorplanetcenter.net/db^earch 
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Three of the new candidates, (311999) 2007 NS2, 2011 SC191 and 2011 
UNg3 have been observed at multiple oppositions, hence their orbits should 
be well determined. The remaining three, 2011 SL25, 2011 SPisg and 2011 
UB256 are only single-opposition discoveries, implying that their orbits are 
uncertain. Therefore we have not considered those further in this paper 
and focus our attention on the three objects with well-determined orbits. 
2007 NS2 was discovered by 0AM Observatory, La Sagra (J75) and was 
subsequently precovered in 1998 LONEOS and LINEAR survey data (MPEG 
2007-009). 2011 SC191 is a recovery, by the Mt Lemmon survey, of 2003 
GX20, an object initially discovered by the NEAT survey. Similarly, 2011 
UNes is a recovery, by the Mt Lemmon survey, of one of its own discoveries, 
2009 SA170. 

Apart from 5261 Eureka, (101429) 1998 VF31 and (121514) 1999 UJ7, the 
study by SMT05 also demonstrated long-term stability for 2001 DH47. This 
object also satisfied the criteria of our search through the MFC database 
but does not appear in the Minor Pl anet Center l ist of M artian Trojans as 



do the other three objects (see also iTodd et al.l . l2012bl ). Interestingly, it 



has not yet been awarded a number despite having been a multi-opposition 
object since 2003. These "anomalies" cast its status as a stable Trojan of 
Mars and/or its orbit in doubt and so it has not been included in Table [1] 
We have nevertheless considered it in the combinatorial analysis (Section [5]) 
insofar as to show that the main conclusions of the paper do not change if it 
is included. 



3. Establishing Trojan stability by long-term numerical integra- 
tions 

To determine (i) whether they are, in fact, Trojans of Mars, (ii) their 
likely lifetime in Trojan libration, and (iii) if not transitory, whether they 
belong to either the FD or SD populations (SMT05), we resort to numer- 
ical integration of their orbits. The covariance ellipsoid for each of these 
three object s and Eureka was p opulated with 101 clones using the method 
described in iDuddy et al.l (120121 ). These clones were then integrated for 10^ 
yr in the past under a force model that included the gravity of the 8 ma- 
jor planets, the dwarf planet Ceres and the large Main Belt asteroid Vesta. 
The code utilises the second-order mixed variable sy mplectic (MVS) algo- 



rithm implemented within the MERCURY package fj Chambers! . 119991 ). A 



model of the along-track component of the diurnal Yarkovsky acceleration 
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(jFarinella et al.l . 119981 ) was also incorporated into MERCURY's user- defined 
force feature. In this model, the magnitude of the Yarkovsky acceleration ay 
has the form 

2 eaT^AT 

^y = —B ^cosC (1) 

pR c T 

where p is the bulk density, R the object radius, e is the surface thermal 
emissivity, a is the Stefan-Boltzmann constant, T the surface temperature 
and c the velocity of light, ( the rotational axis obliquity. AT/T is the 
effective temperature difference around the surface of the body and depends 
on p and R as well as the surface albedo A, specific heat capacity C, thermal 
conductivity K and rotational frequency u. The magnitude of the force 
depends on the heliocentric distance through the term AT/T which is a 
function of solar insolation. In the integrations, the heliocentric distance 
used in Eq. [T]was kept fixed at 1.524 AU, the semimajor axis of the Martian 
orbit. 

Two separate batches of = 101 clones per object were integrated. 
In the first batch the Yarkovsky force was switched off and the system is 
conservative. In the second, each one of N randomly-generated clones were 
assigned a value of the Yarkovsky acceleration ay = ay.max {2i/N — 1), i = 
0, 1, ■ ■ ■ , A^. The value of av.max was calculated from Eq. [Hby setting C = 0, 
an assumed bulk density of 1 g cm~'^ and equal to the surface density, K = 
4 X 10-3 W m-i K-\ C = 680 J kg-^K-\ e = 0.88 and A = 0.12. The 
resulting values oi da/ dt^^^ (= lay^rnax/"^-, n being the mean motion) for the 
four objects are: 80 m yr'^ (5261)^ 150 m yr'^ (311999), 450 m yr'^ (2011 
SCi9i) and 480 m yr"^ (2011 UNgs). 

A rotation period of 6 hr was used for Eureka ( Rivkin et al. . 20031 ) and 



the same value assumed for the other asteroids. For the radius we adopted 



the value of 0.65 km for Eureka (jTrilling et al.l . 120071 ). For the other objects 



the radii were calculated from the absolute magnitude assuming an albedo of 
0.5 using the MFC absolute magnitudes to diameter conversion tablefl This 
is done to err on the side of caution and overestimate the actual Yarkovsky 
acceleration on the asteroids by adopting the smallest size that is consistent 
with their absolute magnitudes. 

In the first instance, we observe that the orbital mean longitude A — Anars 



^ Available at http : / /www.minorplanetcenter .net /iau/lists / Sizes . html 
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and semimajor axis a — ctMars of the new objects relative to Mars librate 
around the location of the L5 equilibrium point (Fig. [T]). The intersection 
of the curves for the smaller amplitude cases is likely due to the signifi- 
cant orbital inclination of Ma rtian Trojans which shifts the libration center 
( iNamouni and Murrayl . |2000| ). 

The results of the full 10*^ yr simulations are presented in Figs [2] and El 
For both batches, the relative longitude A,. = A — AMars (left panels) for all 808 
clones remains within 15° of the L5 equilibrium point. On the right panels 
we have plotted 64-point running means of the osculating eccentricity. In 
the first batch there is little or no dispersion or scatter between the clones' 
end states while significant dispersion is observed for the batch of clones that 
evolved under the Yarkovsky force. The higher scatter for clones of 2011 
UNga and 2011 SC191 in this latter batch is likely due to stronger Yarkovsky 
acceleration (as indicated by the higher values of hmax for these objects) 
although the larger uncertainty in their orbits compared to that of Eureka 
and 2007 NS2 may also play a role. No evidence of gross instability or 
secular growth in the eccentricity is observed. On the basis of this evidence 
we conclude that all three new objects belong to the slow diffusing population 
(SD), able to remain as L5 Mars Trojans for Gyr timescales. 



4. Estimating proper elements of Martian Trojans 

To establish how they relate to each other and the other known Trojans, 
particularly Eureka, we have calculated the proper elements of both new and 
previously known objects (Table |2]) by first carrying out a short (7.3 x 10^ yr) 
integration of their nominal orbits as described above with a 4d time step. 
The integration period was chosen so as to cover at least two periods of the 
z/5 mode (~ 3 X 10^ yr). The output was sampled with an interval of 128d 
to avoid aliasing of Mars' and the Earth's orbital frequencies. The output 
was processed with a programme utilising the Frequen cy Modified Fourier 
Transform algorithm ( ISidlichovsky and Nesvornyl . 119971 ). Our technique for 
estimating the proper amplitu de was somewh at different than that of SMT05. 
It is, essentially, the method of iMilanil ( 1l993l ) with two important differences: 
(i) we use the definitions of d (the proper amplitude of the se mimajor axis) 
and -D (the proper amplitude of the mean longitude) used in iMarzari et al. 
( l2003al Jbl) and SMT05, ie the "full" amplitude from minimum to maximum 
and twic e the v a lue li s ted in proper eleinent ca talogs and corresponding lit- 
erature ( iMilanil . Il993l : iBeauge and Roigi . 1200 ll ). and (ii) the position of the 
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equilibrium point is not fixed at ±60° but is taken directly from the harmonic 
analysis. Our estimate of D for Eureka (~ 11.5°) slightly underestimates the 
value (15°) overplotted in Fig. 3 of SMT05. However, it is in good agree- 
ment with the value used in their Fig. 11. Uncertainties are estimated as 
the differences between the output of this procedure and the results of car- 
rying out the FMFT analysis on the first half of the time series only. As an 
additional test of the robustness of our and /p proper element estimates, 
we calculated 64-point running means of e and I for each clone of the three 
new objects and Eureka at the end of the 100 Myr integrations that included 
the Yarkovsky effect and determine lower and upper bounds after removing 
the five highest and lowest values (10% of the sample). These bracket the 
proper element estimates obtained from harmonic analysis of the short-term 
integrations (Table |2]). In addition, our values of d and estimated inde- 
pendently of each other, satisfy the approximate relationship d ~ y/^JlaD 
where /i and a are the Mars/Su n mass ratio (~ 3 x 10"'') and semimajor axis 
(1.52 AU) of Mars respectively ( Erdi . 1988[ ). The relatively large uncertainty 
for Cp in the case of 1998 VF31 is probably due to its proximity to the 
secular resonance. 

All four objects are contained within the domain 10° < D < 20°, 0.04 < 
ep < 0.075 and 18°. 5 < Ip < 22°. 5 (Table |2]). The corresponding dimensions 
of the stable island found by SMT05 (illustrated in Fig. |4]) can be partitioned 
into 10 such bins in D and 4 bins in each of and I p. The phase space volume 
which contains Eureka and the three new Trojans is then ~ 1/2x10x4x4 = 
100 times smaller than the area of this stable island where the factor of 1/2 
arises due to the roughly triangular projection of the latter in Cp — Ip. 

The new objects were discovered by different NEO surveys and at different 
years and were chosen for this study for the quality of their orbits. Assuming 
a uniform proper element distribution of long-lived Martian Trojans, it is not 
clear why their proper elements should be so near each other and so similar 
to Eureka's. 



Is the clustering of the proper elements of Eureka and the new 
objects statistically significant? 



Cl ustering tests used previously in thi s context (IMilanil . Il993l : iBeauge and Roigl . 



200 ll : IVokrouhlicky and Nesvornyl . 120081 ) were designed to pick out object 
concentrations superimposed on a background population, but such a pop- 
ulation is absent here. To answer the question in our case, we estimate the 
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probability that such a clustering would occur by chance under the assump- 
tion of a uniform distribution of objects in Cp and Ip. The validity of 
this assumption, which may depend on the formation scenario for these as- 
teroids, will be discussed further in the Conclusions section. The projection 
of the domain containing long-lived Trojans (Fig. 1 of SMT05) on the Cp — Ip 
plane may be approximated by a quadrangle (Fig. H]) defined by the points 
(0.02, 13°), (0.02, 28°), (0.15, 28°) and (0.15, 20°). The location of the z/5 
resonance at D ~ 140° indicates that the region of stability may not extend 
as far as Z) = 150°. We proceed as follows: we assume that (a) Martian 
Trojans can occupy a range of at least 100° in D (bottom panel of Fig. Hj) 
and (b) that their distributions in Cp — Ip and in D are uncorrelated. We 
then partition the D domain alone into 10 cells of 10° width each. 

The problem of distributing k consecutively-discovered Trojans into n 
cells is equivalent to seeking the integer, non-negative solutions of the equa- 
tion 



^ ^ Xi — k. 



i=l 



The number of discrete solutions (call it C{n, k, 0)) is 

n + k — ] 
k 

If we impose the "clustering" constraint that Xi > m for at least one i, 
the number of solutions C{n, k,m) is found by assigning Xi = m for one of 
the i s and distributing the remaining k — m units in the remaining n — 1 
cells. Since there are n choices for i we have 

^, J . f n + k - m 
Cin, k.m) = n \ , 

\^ k — m 

and the probability that a cluster of m ( "m-cluster" ) would occur by 
chance is 



P{n, k, m) = C{n, k, m)/C{n, k, 0). 

For n = 10, A; = 5 (all L5 Trojans) and m = 4 we get C(10, 5, 0) = 2002 
and P(10,5,4) = 90/2002 ~ 4.5 x 10"^ or l-in-22 odds of occurring by 
chance. Substituting m = 3 in the above gives 450/2002 ^ 0.22 ie a 3-cluster 
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has 2-in-9 odds of occurring by chance. The algorithm fails for [k/m\ > 1 
(square brackets denote the integer part) since different choices for i no longer 
correspond to foreign sets of solutions. Now we increase the dimensionality 
of the model to include Cp — Ip as well. The above probabilities are reduced 
accordingly e.g. for a 3-cell equal-area partition of the trapezoidal domain 
(Fig. m top panel) the reduction factors are 2/7 for m = 4 and 3/7 for 
m = 3 and the respective probabilities become 1.3 x 10~^ and 9.4 x 10~^. 
Such a partition exists because all four Trojans are contained within the box 
18° < Ip < 23°, 0.04 < ep < 0.075 of surface area < 1/3 of that of the 
trapezoidal domain. The result is robust against adopting more conservative 
bounds for the trapezoidal domain eg bounding Cp at 0.11 instead of 0.15 
and Ip at 26° instead of 28°. 

Note on 2001 DH47 (see Section [2]): According to SMT05 it is an Lg 
Trojan with D ~ 80°. If it is included in our analysis, then k = Q and the 
probability for m = 4 evaluates to P = 450/5005 x 9/28 ~ 2.9 x lO'^. Thus 
the significance of the 4-cluster persists. 



6. Could the clustering of the proper elements of Eureka be due 
to observational bias? 



Observational searches for Martian Trojans are more likely to succeed 
within an area ~ 30° around t he osculating mean lorigitude of each of the two 
triangular equilibrium points (ITabachnik and Evansl . ll999l . l2000bl ). However, 
the new objects were discovered by different NEO surveys and at different 
years, chosen for this study for the quality of their orbits (Section [2]); oscu- 
l ating elements of Troja ns do not relate to proper elements in a linear way 
(IBeauge and Roig| . l200ll ): two out of the three known long-lived Trojans have 
l arge D (> 40°); and a recent search for L5 Trojans yielded no new discoveries 



(iTrilling et al.l . l2006l ). These are strong arguments in favour of (a) the orbital 
clustering of these four asteroids - hereafter referred to as the Eureka cluster 
- being a real feature of the actual distribution of Martian Trojans and (b) 
the discovery of the new long-lived Trojans being the result of the gradually 
increasing sensitivity of small body surveys (ie higher limiting magnitude) 
rather than any observational bias towards Trojans with particular orbital 
properties. 
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7. Comparing the Trojan populations of Mars and Jupiter 



The Jovian Trojan clouds are well populated; the MPC lists 5161 Jupiter 



Trojans as of 21 June, 2012 (http: //www.minorplanetcenter.net/iau/lists/ 



Jupit erTrojans.html) . Their proper elements ( lMilanil . ll993l : lBeauge and Roid . 



200ll ) are distributed throughout the lon g-term stable regions of phase space 



0.15, Ip < 35° and D < 40°) f lMarzari et all . l2003af ). Local con 



centr ations in object dens i ty do exist, but unlike the case of the asteroid 
belt f lZappala and Cellind . Il992l: iMilani and Knezeyid. Il994f) . they are the 
exception rather than the rule (IBroz and Rozehnall . 1201 ll ). To compare the 
Jovian and Martian cases, we partition the Jovian phase space into three- 
dimensional cells of width Acp = 0.05, Alp = 5° and AD = 10°. We then 
read off currently available proper element catalogs one object at a time, 
placing it in the corresponding cell until one of the ten cells in the AD 
partition gains four (4) objects. The number of steps required to achieve 
this (the "score") is then a function of the Cp and Ip coordinates of the 
cell. The values of the proper amplitudes of libration read from the cata- 
logs correspond to D/2 so they have been multiplied by 2 before use. To 
assess the sensitivity of the results to the order that each Trojan is "dis- 
covered", we perform a cyclic shifting of the positions of the first 100 Tro- 
jans (arguably the easiest ones to discover) with a randomly- chosen starting 
point and repeat the process 1000 times. Trials with the first 20 and 50 
Trojans yielded either the same or higher scores. In this work we have con- 
sidered the catalogs of semi-analytically-derived proper elements for 1702 
numbered and multi-opposition Jovian Trojans by C. Beauge and F. Roig 
(jhttp:/ /hamilton. dm. unipi.it/~astdys2/propsynth/petra.pro ) and that of nu- 
merically derived proper elements of 1738 numbered and multi-opposition 
Trojans by Z. Knezevic and A. Milani (superseded in May 2012 by a new ver- 
sion containing 4030 entries; http:/ /hamilton. dm. unipi.it/~astdys2/propsynth/tro.syn). 
The results are almost identical for the two catalogs so we resort to using the 
Beauge and Roig catalog only. We have carried out the procedure in both 
the L4 and L5 clouds and find that the criterion is achieved earlier (ie at a 
lower score) for the L4 cl oud, consistent with the greater density of Trojans 
within the leading cloud (IBeauge and Roigl . I2OOII ) . 

In the left panel of Fig. |5] we show the base 10 logarithm of the score 
achieved for the "nominal" case, ie L4 Trojans are considered in the order 
that they appear in the catalog. The minimum cell score is 52 for the cell 
centered at Cp = 0.075 and Ip = 20°. In the right panel we compare two 
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one-dimensional cross-sections of the data, one for Cp = 0.075 and the other 
for Cp = 0.025 with our result for Mars. Uncertainties are from randomly 
shifting the first 100 Trojan entries in the catalog. The minimum of 36 is 
achieved at the cell mentioned above. Hence, at least 30-40 and typically 
~100 L4 Jovian Trojans need to be "discovered" before one finds a cluster 
with the same spread in proper elements as in the Martian case. From this 
exercise, we conclude that the clustering observed for the first few Martian 
Trojans discovered cannot arise in the Jovian case and that the proper ele- 
ment distributions of the two populations are significantly different. 



8. Discussion: The origin of Martian Trojans 



The existence of a concentration of Trojans in the Martian L5 region 
raises a number of questions regarding their origin and evolution. Cap- 
ture of stable Trojans by the terrestrial planets in the prese nt solar sys- 



tem is unlikely fjSchoU et al.l . l2005l : ISchwarz and Dvorakl . |2012| ) . In one ori 



gin model (A. Morbidelli, reported in SMT05), these objects were captured 
as Trojans of Mars during the planet's random radial wandering in the 
early solar system (cf Fig. 13 of that work). Alternatively, they could have 
been captured into the Trojan regions as a result of planetary mass growth 
( iFleming and Hamiltonl . [2000h . Both scenarios imply deposition early in the 
solar system's history. The three known Trojans are s pectroscopical l y dis 



tinct. Eureka belon gs to the rare A/Sa taxonomic type flRivkin et al.l. 12003 



DeMeo et al.l . 120091 ) not shared by 1998 VF31 and 1999 UJ7 f lRivkin et al. 



20071 ) and suggesting that Eureka was once part of a differentiated parent 



body. Thermal IR spectroscopy with the Spitzer spaceborne facility showed 
a high olivine conte nt, consistent wi th either an oxidised chondritic or achon- 
dritic composition ( Lim et aD . 2011 ). These objects' diverse taxonomic types 
are consistent with radial mixing of the main a steroid belt under a scenario 



requiring the early inward migration of Jupiter (IWalsh et al.l . I2OIII ). If this is 



what actually transpired, the deposition of Martian Trojans predates that of 
Jupiter Trojans, if the latter co incided with the so-called Late Heavy Bom - 



bardment as proposed recently (iGomes et al.l. l2005l: iMorbidelli et al.l. 120051). 



Orbital groupings of asteroids (IZappala and Cellinol. 



1992 



20081 ) can be prod uced either through collisions (iMichel et al 



tational breakout (I Walsh et al 
progenitor body. 



Vokrouhlicky and Nesvornyl . 
20031) 



20081 ) and fission flPravec et all . I2OIOI ) of a 



or ro- 



ll 



CoUisional lifetimes of Main Belt asteroids are a few times 10^ yr for ob 



jects of Eureka's size and a few times 10^ yr for smaller asteroids (jO'Brien 



2004f). Collisions ha ve been important in the evolution of Jovian Trojans 



( Marzari et al.l . ll997l ): Trojan- Trojan collision probabilities are higher (jdeirOro et al 



19981 ) than those in the main belt due to the spatial confining action of the 
resonant libration. No studies directly relevant to Martian Trojan coUisional 
evolution exist to-date but the two most important impactor populations will 
likely be (a) Mars Crossers (b) other Martian Trojans. In the latter case, the 
efficiency of coUisional comminution would depend on the (as yet unknown) 
population of small Trojan impactors. Due to the high inclination (~ 20°) 
of the Trojans and their proximity to the Sun, typical collision speeds with 
both populations would be more sim ilar to those between MBAs and NEAs 

> 



or NEAs and NEAs, > 10 km s"^ fiBottke et al. 



19941 ). 
1.524 AU, e 



For example, at 
= 0, / = 20° and 



the intersection point of two orbits with a 
AQ = 90°, the relative speed is ~ 11 km sec~^. Collision probabilities may 
be enhanced relative to the Jovian case as the spatial azimuthal extent of 
the Trojan clouds scales as a~^, the semimajor axis 



Studies of the coUisional fragmentat i on of asteroids (lAsphaug et al.l . 11999 ; 
Benz and Asphaud . Il999l : iMichel et al.l . |2003| ) apply mainly to lo wer velocity 



impa cts and larger target bodies. In recent hydrocode simulations (jJutzi et al. 



2010l ). production of large fragments is enhanced at high impact velocities or 
km-sized targets while ejection velocities are a few m s~^ at most for frag- 
ments > 1/10 the size of the target. The resulting semimajor axis change 
of a fragment must be smaller than the radial width of the Trojan region 
(~ ay/8i j/3 where /x is the planet's rn ass parameter and a its orbital semima- 
jor axis flMurray and Dermottl . Il999l )). The corresponding change in along- 
track velocity magnitude for an object in a circular, planar orbit about the 
Sun is 



Av = -riA/S/i/Sc 



(2) 



where n is the planet's orbital mean motion. Equation |5]evaluates to ~ 10 
m s~^ for Mars and ~ 350 m s~^ for Jupiter. Fragment size is inversely pro- 
portional to ejec tion velocity and largest fragment ejection velocities increase 
with target size (jJutzi et al.l.l2010[). This has tw o important consequences: (i) 
unlike the case of Jupiter ( iMarzari et al.l . 119971 ). all but the largest coUisional 
fragments of Martian Trojan progenitors would have sufficiently high veloc- 
ity to directly escape the Trojan regions, resulting in the observed paucity 
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of Eureka-sized objects, and (ii) the Eureka cluster's Trojan precursor could 
not have been larger than a few tens of km acros s, such as th o se po stulated 
in some models of early giant planet migration f Walsh et al. . 201ll ). For a 
mass ratio between parent body and largest fragment of ~_0.5 (i.e. near 
the catastrophic disruption threshold (IBenz and Asphaugi . Il999l )) the num- 
ber of collisions required to produce an object of Eureka's size (r = 600 m) is 
10 log R + 2.2 where R is the parent body's radius in km. Setting i? = 5 — 50 
km yields 9-19 collisions over the age of the solar system, implying 4-8 colli- 
sions over the past 2 Gyr that the cluster likely formed (see below). We note 
tha t NEA-NEA c ollisio nal lifetimes for objects 1-30 km across are < 10^ 



yr (IBottke et al.l . 119941 ) but those in the Martian Trojan clouds could be 



significantly shorter a s discussed above. 

Rotational fission (jjacobson and Scheeresl . 1201 if ) has been cla i med t o be 
responsible for the so-called "spin barrier" for NEAs (iPravec et al.l. 120071) and 



pairs of Main Belt asteroi ds with unusually similar orbits fIVokrouhlicky and Nesvorny 



20081 : iPravec et al.l . l2010l ). Characteristics of the mechanism are a separation 



velocity comparable to the escape velocity, ~ 1 m s~^ for a 1 km progenitor, 
and mass ratios < 0.2, implying a size ratio < 0.6. The latter constraint is 
satisfied for all possible pairs among these four asteroids that include Eu- 
reka and pairs that include 2007 NS2 and the two smaller objects. Sur- 
face escape velocity can be expressed as a function of the object's radius r: 
Vesc = ^/S/STrGpr"^ where G is the gravitational constant (~ 6.67 x 10~^^ N 
m^ kg~^) and p the bulk density. If the seco ndary fissions away from the pri- 
mary at escape velocity (IPravec et al.U2010l ). r must be < 8 km for p = 3000 
kg m~^ to prevent escape from the Trojan region (see previous paragraph). 
The fundamental difference between this mechanism and the collisional one 
is that ejection speed is not a strong function of ejectum size; therefore, even 
small secondaries (mass ratio -C 0.2) should survive as Trojans for a small 
enoug h primary. The only available estimate of Eureka' s rotation per i od (> 
4.8 hr: iRivkin et al.l . l2003l ) is not near the "spin barrier" (IPravec et al.l . |2007| ) 
but is consistent with the expected prima r y rota tion rate for the mass ratios 
of revelance here (jJacobson and Scheeresl . I2OIII ) . We note that phase cover- 
age of Eureka's lightcurve is incomplete; additional photometry is desirable 
to confirm the rotation rate value. 

Both mechanisms predict low (< 10 m s~^) separation velocities for the 
three new Trojans if they originated from Eureka or a prog e nitor of compa- 
rable size. Applying the orbit comparison metric by iMilanil (Il993[ ). given by 
the expression 
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na 



y^(l/4) {6d/af + 2 {6epf + 2 {6 si 



sin /„ 



(3) 



we find that the magnitudes of mutual orbital differences vary from 350 
m s-i (2007 NS2 - 2011 UNgg) to 1800 m s"^ (5261 - 2011 SC191). Hence, 
some orbital diffusion must have taken place within the Eureka cluster to 
modify Cp and ip. This can be used to constrain the time of formation. 

The proper elements of long-lived dynamical clones of Eureka diffuse ap- 
preciably over the age of the solar system (SMT05). The proper eccentricity 
diffuses faster, reaching values up to 0.15 over 2 Gyr. During the same period, 
the proper libration amplitudes of the clones remain similar to, and slightly 
larger than, Eureka's. They spread over a range of ~ 50° over the next 2 Gyr. 
The objects investigated here have 0.04 < Cp < 0.075 and 10° < D < 20° and 
their libration amplitudes are slightly larger than Eureka's. Consequently, if 
these orbital differences are due to the action of chaotic diffusion alone (but 
see paragraph below), the event or events that produced them have taken 
place less than 2 Gyr ago. This is consistent with a formation relatively late 
in the precursor's coUisional history (ie a parent body originally tens of km 
across but close to Eureka's present size when the cluster was formed) in 
order to keep ejection velocities low. 

An additional complicating factor may be the action of the Yarkovsky 



Nesvornv et al.. 


2002 


) . It is thought 


modification ( 


Moldovan et al.. 


2010; 



Broz and Rozehnall . I2OIII ) but may be more important in Mars' case since 



for reasonable choices of asteroid physical parameters (see lBottke et al. 



20061 . and references therein). Dissipa tive forces in g eneral modify the linear 
stability properties of Trojan motion (jMurrayl . 119941 ) in a way unique to the 
force considered. SMT05 found that the Yarkovsky effect depletes the Mar- 
tian Trojan regions of objects < 5 m in size over the age of the solar system. 
In this work, we found that the Yarkovsky force affects the long-term evo- 
lution of the eccentricity. That such a size- and rotational state-dependent 
force has not dispersed the Eureka cluster beyond the point of recognition 
implies that either (a) it it inefficient, at least on the timescales of the clus- 
ter's age, or (b) it is similarly efficient on the members of the cluster. The 
latter possibility is unlikely given the disparate sizes of its members and their 
presumably different (but unknown) rotational axis orientations. We cannot, 
however, exclude the possibility that the high relative velocities (> 300 m 



14 



s~^) between the orbits of cluster members are due, at least partly, to this 
effect. In that case, the dynamical age of the cluster would need to be re- 
vised downwards with respect to models dealing with chaotic diffusion alone. 
Although beyond the scope of the present study, follow-on work to quantify 
the long-term effects of the Yarkovsky force on Martian Trojans would be 
desirable. 

An additional formation model is worthy of consid eration here, that of the 
direct coUisional capture of Troj ans as proposed by [Colombo and Franklin 



( I1971I ) and recently explored by iTurrini and Marzaril ( l2009l ) in the context 



of the origin of the irregular satellites of Saturn. In this scenario, a collision 
between two passing asteroids generates a spray of fragments and, among 
them, some have the right velocity to become Trojans. To explain the ob- 
servations, such an event would have to produce objects within the limited 
phase space volume occupied by Martian Trojans. This implies a fragment 
velocity field such that, in the first instance, the heliocentric semimajor axis 
dispersion is no more than a few times 10^^ AU. Indeed, it is not clear why 
such an event would not fill the stable island identified by SMT05 with frag- 
ments but instead produce the concentration of Martian Trojans reported in 
this paper. In conclusion, although such a scenario is possible, its plausibility 
remains to be demonstrated. 



9. Conclusions and Implications 

This work has demonstrated the following: 

1. The existence of additional long-lived Trojans of Mars 

2. A statistically significant concentration of Trojans around the largest 
known object of this class, 5261 Eureka 

3. A significant difference between the orbital distributions of Trojans of 
Mars and Jupiter 

A potential caveat on conclusions 2 and 3 is our assumption of an ini- 
tially uniform distribution of Trojans. However, if the capture mechanism in 
the case of Mars produced an initial non-uniform distribution within the is- 
land of stability discussed in Section |U the action of chaotic diffusion would 
have displaced these objects from their original locations and diluted any 
primordial features of their distribution over the age of the solar system. Al- 
ternatively, one could argue that the rate of diffusion varies considerably for 
different position within this stable island but the results of SMT05, showing 
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similar lifetimes for different long-lived Martian Trojans, do not support such 
a hypothesis. 

Consequently, we have examined the following hypotheses for the exis- 
tence of this cluster: 

1. That it is the result of the collisional fragmentation of a progenitor 
body. This scenario is consistent with current understanding of asteroid 
collisional evolution and current knowledge of these Trojans' physical 
properties. No studies specifically applicable to Martian Trojans exist 
to-date. 

2. That it is the result of the rotational fission of a progenitor body. This 
is also consistent with our understanding of the process and available 
physical information. In addition, it constrains the - largely unknown 
- rotational states of the asteroids. 

3. That it is the result of direct collisional capture of asteroids in the 
Trojan region. We consider this a possible, but not yet a plausible, 
scenario pending a quantitative demonstration of its efficiency. 

For the first two scenaria, the mutual orbital dispersion between members 
of this Trojan cluster implies an upper age limit of 2 Gyr. 

New observational surveys will soon extend th e known population of Mar- 
tian and Terrestrial Trojans to < 100 m objects fiTodd et all l2012al lbh. The 



orbital distribution of new discoveries will constitute a test of the relative 
importance of the processes discussed above. If collisions are dominant, the 
paucity of km-sized Martian Trojans should be at least as apparent, if not 
more, at the smaller sizes. On the other hand, rotational fission will produce 
small asteroids that will remain within the Trojan clouds, shoring up the 
size distribution of the population at the small end. In effect, Trojan dy- 
namics become a "sieve" to distinguish between different models of Martian 
Trojan production. Due t o the unique geometri c const raints of observing 



Terrestrial planet Trojans (ITabachnik and Evansl . l2000bl ). the case of Mars 



is the most suitable to serve this function. In any case, interpretation of the 
distribution of the newly discovered Martian Trojans, especially at the small 
end of the population, should be accompanied by (a) a re-examination of 
non-gravitational forces under the unique circumstances of Martian Trojan 
dynamics and (b) observational characterisation of key physical properties 
- size, albedo, taxonomic type and rotational state - of known objects and 
future discoveries. Finally, searches for genetic asteroid groups in the Terres- 
trial planet region are compromised by frequent planetary close approaches 
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(IFu et al.l . I2OO5I : I Jacobson and Scheeresl . l2011uSchunova et al.U2012l ). These, 



by definition, do not affect stable Trojans of the Terrestrial planets. In this 
sense, our work highlights the Trojan regions of the inner solar system and, 
in particular, the Martian Trojan clouds as unique "isolation chambers" to 
identify asteroid production and family creation events close to the Earth 
and gain new insight on the physical and dynamical evolution of our solar 
system's small bodies. 
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List of Figures 



Libration around the corresponding equilibrium point (plus 
sign) for all six long-lived Martian Trojans considered in this 
work. In order of decreasing libration amplitude, they are: 
(121519) 1999 UJ7 (L4; dashed blue curve), (101514) 1998 
VF31 (L5; blue curve), (5261) Eureka (L5; red curve); (301999) 
2007 NS2, 2011 SC191 and 2011 UNgs (L5; amber curves). . . . |28 
The dynamical evolution of 404 clones (101 clones per object) 
in our 10^ yr integrations of the fully conservative system. Left 
panels: Longitude relative to Mars A — Anars vs relative nor- 
malised semimajor axis (a — CMars) /^Mars- Only one clone in 
every ten is shown for clarity. Right panels: 64-point running 
means of the osculating eccentricity vs time. Only one clone 



in five is shown for consistency with Fig |3l |29 

As for Fig. [2] but including the model of the Yarkovsky force 
described in Section [31 Left: Only one clone in ten is shown 



for clarity. Right: Only one clone in five is shown for clarity. . |30 
Domains in proper element space shown by SMT05 to host 
the most stable Martian Trojans. These are enclosed within 
the dashed lines. Grey curves indicate the locations of secular 
resonances discussed in that work. The coloured ellipses in- 
dicate the location of the six Trojans and their semiaxes the 
uncertainties in their proper element estimates. Blue dashed 
curve: (121514) 1999 UJ7, Blue curve: (101429) 1998 VF31, 
Amber curve: (5261) Eureka, Brown curves: (301999) 2007 
NS2, 2OIISC191 and 2011UN63. Top panel: Proper eccentric- 
ity (cp) vs proper inclination (Jp). For secular resonance plot- 
ting, the proper amplitude of libration D was fixed at 40°. 
Bottom panel: Proper amplitude of libration vs proper eccen- 
tricity. For secular resonance plotting, Ip was fixed at Ip = 20°. 
(301999) 2007 NS2 and 2011 UNgs have been moved slightly 



farther apart for clarity. |31 
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Left panel: Logarithmic density plot of the score achieved, by 
each resolution element in proper element space, of the popu- 
lation of known Jovian Trojans under the procedure described 
in the text. Each resolution element corresponds to objects in 
the 2-D domain (cp - 0.025, Cp + 0.025) x (Jp - 2°.5, Jp + 2°. 5). 
The open black circle indicates the location of the Mars Tro- 
jan cluster identified in this work. Right panel: Section of 
left panel for ep = 0.025 (filled grey boxes) and Cp = 0.075 
(filled grey circles). The red filled circle represents the result 
obtained in this work for Mars 
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Tabic 1: Orbital elements and physical properties of asteroids considered in this work. 



Designation 


Epoch 


a 




T 
I 




LJ 






(JD) 


(AU) 


e 


(deg) 


H 


(km) 


Opps 


(5261) Eureka 


2456000.5 


1.523482 


0.06483 


20.28 


16.01 


1.28 


13 


(311999) 2007 NS2 


5) 


1.523705 


0.05399 


18.62 


17.75 


[0.67-2.1] 


3 


2011 UN63 


)) 


1.523709 


0.06471 


20.36 


19.73 


[0.22-0.66] 


2 


2011 SC191 


5) 


1.523811 


0.04406 


18.75 


19.44 


[0.25-0.77] 


3 


2011 SL25 


2455800.5 


1.523888 


0.11455 


21.50 


19.5 


[0.24-0.74] 


1 


2011 SPi89 


2455840.5 


1.523279 


0.04040 


19.88 


21.2 


[0.11-0.33] 


1 


2011 UB256 


2455860.5 


1.524049 


0.07118 


24.31 


20.1 


[0.18-0.57] 


1 


(101429) 1998 VF31 


2456000.5 


1.524213 


0.01002 


31.30 


17.0 


0.78 


6 


(121514) 1999 UJ7 


2455800.5 


1.524445 


0.03915 


16.74 


16.9 


[0.78-2.3] 


6 



The diameters of 5261 and 1998 VF31 are from lTriUing et aLl ( 120071 ). For the other objects we 
used the MFC's Absolute Magnitude to Diameter conversion tables to produce likely diameter 
ranges for albedos between 0.05 and 0.5. 
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Tabic 2: Proper elements for all long-lived Martian Trojans from our simulations. 







d 


D 








< / > 


Designation Region (xlO^AU) 


n 




n 


< e > 


n 


(5261) Eureka 


L5 


2,34±o.o2 


11.4±o-i 


0.0570°2i 


22.20=^°-25 


0.060°5 


22.35^0-25 


(311999) 2007 NS2 


L5 


3_54±o.oi 


13.9±o-i 


0.0444°°3 


20.96±°-i^ 


0.045°^ 




2011 UNea 


L5 


34Q±o.o9 


13.8±°-^ 


0.0466°2i 


21.58±°-28 


0.0525^25 


21.6±o-9 


2011 SC191 


L5 


4_42±o.i5 


17_4±o.5 


0.0717°2i 


19.15±o-52 


0.071°6 


19.15^°-65 


(101429) 1998 VF31 


L5 


g_74±o.80 


38.1±2-4 


0.0817^15 


32.o7±oo3 






(121514) 1999 UJ7 


u 


17_30±o.o9 


70.7±o.3 


0.0346°^^ 


18.05^°-6i 






Columns 3, 4, 5 & 6 


give 


the proper amplitudes of the semimajor axis 


and mean 


longitude 



with respect to Mars, proper eccentricity and proper inclination from our FMFT analysis. 
The 7th and 8th columns give the range of values spanned by 64-point running means at the 
end of the 10^ yr integrations as explained in Section [31 Uncertainties are the superscripted 
values except for columns 5 and 7 where they replace the bolded digits. 
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